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Abstract 

One may impose to a system with spontaneous broken symmetry, boundary condi- 
tions which correspond to different pure states at two ends of a sample. For a discrete 
Ising-hke broken symmetry, boundary conditions with opposite spins in two parallel 
limiting planes, generate an interface and a cost in free energy per unit area of the 
interface. For continuum symmetries the order parameter interpolates smoothly be- 
tween the end planes carrying two different directions of the order parameter. The 
cost in free energy is then proportional to L"'"^ for a system of characteristic size L. 
The power of L is related to the lower critical dimension, and the coefficient of this 
additional free energy vanishes at the critical temperature. In this note it is shown 
within a loop expansion that one does find the expected behavior of this twist free 
energy. This is a preamble to the study of situations where the broken continuum 
symmetry is believed to be more complex, as in Parisi's ansatz for the Edwards- 
Anderson spin glass. 
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1 Introduction 



Spontaneously broken symmetries are characterized by the existence of several pos- 
sible pure states. If one imposes "twisted" boundary conditions , i.e. different pure 
states at two ends of the system, the free energy per unit volume will be slightly 
greater than the free energy corresponding to one single pure state over the whole 
system. 

For a simple discrete symmetry, such as the Z2-symmetry of Ising-like systems, one 
may consider an (hyper)-cubic system with up spins in the z = plane, down spins 
in the z = L plane and for instance periodic boundary conditions in the transverse 
directions Xi,X2,- ■ - xa-i- This will generate an interface in the system centered 
around some plane z = Zq and a cost in free energy 

AF = - = aL"-' (1) 

in which the interfacial tension cr(T) is finite at low temperature, but vanishes at the 
critical temperature as 

cr{T) = ao(^^)\ (2) 

As is well-known the power (d—l) of L in ([^) implies that the lower critical dimension 
of systems with a discrete symmetry is equal to one, i.e. there is no ordered phase 
unless d is greater than one. Widom has first proposed the scaling law 

X = {d-l)u (3) 

in which u is the correction length exponent 

e = eo(^^r- (4) 

The corresponding amplitude relation implies that the combination 

^o'-Vo (5) 

is universal. All this was studied long ago by renormalization group techniques 
and (4 — c?)-expansion. At leading order the classical (mean field) solution is a kink. 



of hyperbolic tangent shape, interpolating between up and down spins, and fluctua- 
tions are given at one-loop order by the Fredholm determinant of a one- dimensional 
Schrodinger operator in a l/cosh^(z — zq) potential which, as is well-known, is solv- 
able analytically. 

For more complex spontaneously broken symmetries, continuum symmetries, or 
replica-symmetry breaking, the situation is less trivial, and it it is necessary to look 
into the problem in some detail in order to understand the lower critical dimension. 
For a continuum symmetry group G, broken down to a subgroup H, as in N- vector 
models, one considers the free energy with two different pure states in the two planes 
z = and z = L. For an N- vector model one considers for instance an order parameter 
uniform along the vector (1, 0, ■ ■ ■ , 0) in the z = plane, and uniform but rotated 
by an angle 6*0 in the plane z = L, i.e. lying along the vector (cos 9o, sin 6'o, 0, ■ ■ ■ , 0). 
There again one expects a cost in free energy 

AF = a{T, 60)1'^-'^ (6) 

in agreement with a lower critical dimension equal to two, and with a "twist" energy 
cr{T,6o) (or spin stiffness constant) vanishing as 6^ for small ^0;( the ratio (j/6q is 
the helicity modulus defined by Fisher, Barber and Jasnow 0), and vanishing at 
0, ^ like (Tc — Ty^'^~'^\ If it is quite elementary to verify these statements within 
mean field theory, not difficult also to check them in the vicinity of the lower critical 
dimension di = 2 through the non-linear sigma model |^, it is not so simple to 
examine the problem below the upper critical dimension du = 4. This note is thus 
devoted to this point. Our aim in performing this calculation is to repeat it later 
for a spin glass . There, below the temperature of transition, one recovers a broken 
continuum symmetry, displayed by Parisi's ansatz P] of replica symmetry breaking 
for the Edwards-Anderson model, which yields a continuum of schemes with equal 
free energy (reparametrization invariance). If the situation at du = Q is more or less 
under control, despite its complexity , the knowledge about di is poor, although it 
is believed to lie between two and three |T^, . If one imposes again two different 
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schemes at two ends of the system, one expects W^ a cost in free energy 



AF = (7) 

with some negative anomaly r] which would yield a lower critical dimension di = 
2 — 1]. The possible presence of this anomaly requires to compute at least a one- 
loop correction to mean field theory. This difficult calculation will be reported in a 
subsequent article, and this note simply aims at showing that already for the well 
understood N-vector model, the theory is somewhat involved. The rest of this note 
is thus devoted to the N-vector model below four dimensions, treated thus through 
a (0^)^ field theory and an e = A — d expansion.lt is interesting to note that a 
direct calculation of the helicity modulus has also been performed directly in three 
dimensions, in spite of the singularities expected from Goldstone massless modes 
1151 , p!6[| . This has been done by keeping a symmetry breaking field until one can let 
it go safely to zero at the end of the calculation of this helicity modulus. 

2 Mean field theory 

The action for the N-vector model in the broken symmetry domain is 

S = J\z J d'-W [^(V0)^ - ^\m' + (8) 

in which t is proportional to T — Tg. A pure state throughout the bulk would have a 
magnetization M whose magnitude is given by 

\t\ = gM\ (9) 

Subtracting the bulk contribution one thus has 

= j^dz J d'^-'x^^iV^f + - M^]. (10) 

The free energy AF is the value of the minimum of A^* with the boundary conditions 

0(z = O,x^) = M(l,0,---,0) 

(f){z = L,x±) = M(cos^o,sin^o,0, ■■■,0). (11) 



We fix here the value of the order parameter on the edges, rather than imposing 
magnetic fields on the boundaries. Our partition function will thus be defined with 
fixed prescribed values of the order parameter on the two edges, rather than fixing 
a surface magnetic field and letting the surface order order parameter fluctuate, as 



in the work of M. Krech [T^ for instance. For an = 1 (scalar) order parameter, 
we would have to fix the surface order parameter to a value slightly smaller than the 
bulk magnetization, but for > 1 one can directly take the modulus of the surface 
order parameter equal to the bulk magnetization, as shown in the mean field solution 
of the equations of motion below. It is easy to verify that AS is minimum 

• when the order parameter remains in the 2-plane of the two vectors 
defined by the boundary conditions 

• when is a function of z-alone, i.e. independent of x±. 
and one can parametrize the mean field solution as 

(p = p{z){cos9{z),sm9{z),0, - ■ ■ ,0), (12) 

for which 



The solution will be close to that of an order parameter uniformly rotating between 
the two planes with a constant magnitude M, namely p{z) = M and 6{z) = ^^o, for 
which AF = ^6q^M'^L'^~'^. However, although the solution is close to that for large 
L, we shall need the corrections of order 1/L^ to that simple ansatz, and one has to 
solve the variational equations 



dz 



Up'o') = 



p" - pe'^ - gp{p' ^ M') = 0. (14) 
Defining the dimensionless variables 



1 
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the equations of motion are easily cast into the form 

. dr 

T 

,de 



with 



r^— = V7, (16) 



vir) = ^^-\il-rr- (17) 



We can think of the equation for r as an equation of motion in (r, r)-plane in which 
r starts at r = 1 for r = , decreases down to some Tq , then increases and returns 
to r = 1 at 



To = L^2gM^ . (18) 

The parameters 7 and u) have still to be determined by the boundary conditions. Of 
course the exact solution of the equations of motion (|1^) involves elliptic functions 



13]. However it turns out that it is sufficient for our purpose to consider the regime 



in which 7 is small, which corresponds to L large compared to the correlation length 
^ (or if L/^ is finite, corresponds to small 6q). Indeed in that regime the order 
parameter has essentially a fixed length r(r) close to 1, and 



2gM'^L^ 

uj = 7 + 0(7^) 

r = 1 -7s(r) +0(7^). (19) 

The full integration to this order in ^ — -is then easy and leads to 

, , cosh 1^ — r| , , 

To that same order one finds 

AF = -QlM'^L'^-'' + 0(7^) = —QlL^~''\t\ + 0(7"). (21) 
2 2g 
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Let us note that, in mean field, the correlation length is related to the temperature 

by 

= 2gM'^ = 2\tl (22) 

and thus 7 is small either because L/^is large or because 9q is small. The result (^T]) 
is thus in agreement with our expectations 

AF = aL'^-'^ (23) 

with a = ^^oW vanishing at the critical temperature . We also verify the scaling 
law a{t) ~ (which is expected to be true for d < 4) in four dimensions at 

which z/ = 1/2 and z/((i — 2) = 1. 

As far as mean field theory is concerned the picture is simple : for 6qC,/L small, 
the magnitude of the order parameter remains close to M over the whole sample, 
and its direction smoothly interpolates between the two end planes with a constant 
angle gradient. If we went beyond this simple picture in (^0]) it is because this will be 
needed in the loop expansion when we consider fluctuations around the mean field. 

3 One loop corrections 

We now go to dimension d = A — e and work to first order in e, which requires the 
calculation of one-loop fluctuations around mean field theory. Instead of an ultra- 
violet cut-off given by some lattice spacing, it turns out to be much more convenient, 
as often, to use dimensional regularization. The mean field solution is 

0c = Mr{z){cos 9{z), sin e{z), 0, ■ ■ ■ , 0), (24) 

with r and 6 described in the previous section. It is convenient to introduce an 
orthonormal moving frame consisting of the vectors 

6*1 = (cos^(z), sin^(2;), 0, ■ ■ ■ , 0) 

62 = (-sin^(z),cos^(z),0,---,0), (25) 
6 



plus the {N — 2) fixed unit vectors e^, (a = 3, • • • , N) perpendicular to the two-plane 

— * 

(1-2). The field (l){z,x±_) is then parametrized as 

N 

4){z, afL) = {p{z) + i^iiz, x1))ei{z) + ip2{z, xl)e2{z) + ^ ipaiz, x±)ea. (26) 

a=3 

The boundary conditions on those ipa are periodic in the tranverse directions and, 
since the mean field order parameter phic is equal to the magnetization on the bound- 
aries, one has to impose Dirichlet conditions on the fluctuating flelds ipa{z = 0) = 
ipa{z — L) = 0. A one-loop calculation requires to keep only the quadratic terms in 
■0a of the action. Collecting those terms one finds 

S^So + S2 (27) 

in which 5*0 is the mean field contribution and 

S2 = j^dzj d'^-W[^ EWaY + ^i2gM' + + 3gM\r\z) - iM 



Tz^"^'^ - ^^^^ + -2^^d-z^ + - '^^^^ 



1 ^ 

+ -^M^(r^(^)-1)EV'^]- (28) 

^ 3 

The one-loop free energy is thus equal to the properly normalized 



AF^^o + ^IVln - (29) 

The normalization will be chosen such that AF vanishes with 6q. 
Contribution of the N fluctuating modes 

• The transverse modes '0a, a = 3, • • • , are decoupled and give a con- 
tribution to AF equal to 

]-{N - 2)Trln[-V^ + gM^{i^{7) - 1)] 

2)1''-^ J ^^Tr ln[ql - ^ + gM2(r^(z) - 1)](30) 



7 



On should note that although r'^{z) — 1 is negative, the spectrum of 
— ^ is bounded below by vr^/L^ since we have Dirichlet boundary 
conditions on the planes z = and z = L. Taking the explicit solution 
( p!9| , pOD one sees that the spectrum of —-^ + gM'^{r'^{z) — 1) is bounded 
below by (tt^ — O'D/LP' and is thus positive. 

Therefore a priori one has to compute the Fredholm determinant of 
a one-dimensional Schrodinger operator in the complicated potential 
r'^{z). However for large L, perturbation theory gives very simply the 
answer since r'^{z) — 1 is of order This is to be contrasted with a 

localized Ising interface, for which there is no small parameter for large 
L. The simplification here is due to the fact that the order parameter 
turns slowly from one end of the system to the other one and thus has 
only small fluctuations in the moving frame that we have introduced. 

Then we may replace Trln[q5_ ~ + gM^(r^(z) — 1)] (subtracted to 
vanish at 9q = 0) by gM'^TT[q\ — J^]~^(r^(z) — 1). Expanding on the 
basis of the Dirichlet eigenstates of — the states -y/^sin (nvrz/L), 
we obtain the contribution of these modes, in the large L limit, under 
the form 



in which s{z) is the explicit mean field correction (pOl) . In the large 
L limit one can replace s{z) by one, the sum over n by an integral 
which, combined with the integral over gives the integral L J d'^p/p^ 
which vanishes in dimensional regularization. Those modes have thus 
a vanishing contribution to the terms proportional to L'^"^ of AF. 

• We now come to the coupled 'ipi-ip2 modes, using again that r^(z)— 1 and 
{d9/dzY are of order This allows one again to use a perturbation 

expansion about a massless ?/'2-niode and a massive ■^i. After a lengthy. 
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but elementary calculation, we obtain the contribution of these two 
modes to AF under the form of a sum of five terms : 

o2 1 1 

^^one-loop _ rj^/ |_ 



2 • 



2L2 '-V2 + 2gM2 -V 

%Tr(3--^-— s(z) + 4-s(z)) 



2L2 ' -V2 + 2gM2 ' ' -V2 

^2 1 

2L2 ^(-V2 + 2gM2)(-V2)^ dz^ 



We leave the detail of the calculations to an appendix and simply report 
the result. We have computed the 1/e pole of this expression, for 
arbitrary L/C, and obtained. 

AF = + ^^elL''~'\t\'-'/^ (32) 

(we have kept it under this form since g and have the same di- 
mension). In this expression we have kept the finite (|t| In |t|) term and 
neglected the non-logarithmic terms. 

4 Renormalization and scaling 



We first note that the pole in 1/e in (|32D is independent of L/C, as it should, since 
the renormalizations are independent of this ratio. Next we note that the limit 
of e going to zero should be finite, provided we perform a coupling constant and 
mass renormalization (there is no wave function renormalization at this one-loop 
order). Taking the standard one loop result from literature |]T^ ( with the appropriate 
normlization of the coupling constant chosen in (H)) one has, at one-loop. 



1 1 N + 8 

- = f^ ( 

g 9r Svr^e 



/^-^(- - ^) (33) 



for the coupling constant renormalization (/x is an arbitrary inverse length scale) and 

t = tn{l + gJ^) (34) 
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for the mass (i.e. temperature) renormalization. This gives a renormahzed expression 
for AF in terms of and tR which is finite, as expected, when e goes to zero: 

A^.|-i-Vl<.|(l-||ln^). (35) 

The scahng of the coefficient of 9qL'^~'^ in the critical region requires a replacement 
of qr by the infra-red stable fixed point 

9R-^ + 0{e% (36) 



and the exponentiation 



(l-^lnM)^(M)3.y8.^. (37) 



Given that the correlation legth exponent v has the expansion 



one verifies to this order that 

l-|| = ^(rf-2) (39) 

which does yield the expected scahng law for the vanishing of the twist energy at Tc 
in the 0(A/")-model. 

5 Final remarks 

Although a priori more cumbersome than the calculation of the interfacial energy 
for a discrete symmetry, it turns out that, for a continuum symmetry, it is possible 
to compute the complicated Fredholm determinant of fluctuations around mean field 
theory by an expansion in powers of 1/L which was not available for an interfacial 
wall. The calculation involves a description of the mean field solution in which it 
is not sufficient to simply assume that the order parameter rotates with a constant 
angle gradient from end to end, with a fixed length equal to the magnetization. The 
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calculation presented here may be easily generalized to any continuum symmetry 
group G, broken down spontaneously below a critical temperature to a subgroup H, 
with an order parameter in a given irreducible representation R of G. 

This calculation provides an explicit test of the fact that the renormalizations 
are the same around any background solution: in the usual case one expands about 
a classical solution which is constant over the sample, whereas here one expanded 
around a non-trivial solution, and yet we found that the same coupling constant 
and mass renormalizations did work. We have also verified that the finiteness of the 
end result of the free energy for any ratio L/^. However the method that we have 
followed, has made use of a small parameter, namely 9o^/L. Away from the critical 
temperature this parameter is small because L is large. However if L/^ is finite our 
calculation is restricted to small ^o- In the finite L/^ regime, AF is a priori a more 
complicated function of 60 for which we have only determined the first term. Let 
us stress also that we have used the e-expansion, since we wanted to determine the 
behavior of the twist free energy near the upper critical dimension. 



11 



Appendix : One-loop divergences 



Let us return to the five terms contained in (^Tj) for the one-loop calculation of 
AF. 



fl2 1 

— Tr( ). (1) 



' ' 2L2 '-V2 + 2gM2' 
If L goes to infinity first we may simply neglect the quantization of the 
longitudinal modes and write 

in which it is understood that we have neglected the terms of order 
e°. For finite L/^ the calculation is much more involved. Going to 
the large L limit for the tranverse periodic directions, but keeping the 
quantization of the longitudinal modes one has 



^L''-\2gM^Y-''^ r dx f'"' [coth/y/^^- / ], 

(3) 



in which 

l = L^2^^ = L/i. (4) 

We now use the identity 

poo x'^~'^ 1 

/ dx , [coth/Vx^ + l 1 

/■oo 2 1 

= / [1 _ _1 

/•oo x'^~^ 1 1 

+ / dx , [coth/Vx^ + l + — 1. (5) 
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The first term of tlie r.li.s. of (^) is elementary and gives — l/(2e) plus 
finite terms. It is easy to see that the second integral of the r.h.s. of 
@ is finite when d A. This proves that the divergent part of (a) is 
as expected independent of I = L/^. 



Again if L goes to infinity first 

which vanishes in the dimensional regularization scheme. 



= -2Z^^^(-V^ + 2gM2^(^)) 

3gg L^-^ /• 2 ^ sin^nTTz/L 

2L2 (27r)'^-i y ^^Lio Y^i + ^'^V^' + 2^/M2^ 

Again if one lets L go to infinity first one can replace s{z) by one, the 
calculations are then elementary and yield 

(c) = -^{2gM'y-^/'elL'^-\ (9) 

For finite L/C, one can prove with the help of the explicit form for s{z) 
that the divergent part is unchanged. 



{d) = -^Tr(-L_s(z)) (10) 
Again it is easy with the same integral representation 
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to prove that the leading term, proportional to L'' ^, multiplies the 
integral / d'^p/p^ which vanishes. Therefore 



{d) = (12) 



fl2 1 pi2 

^' 2L2 ^(-V2 + 2gM2)(-V2)^ dz^'' 

(13) 

which, in the large L limit, goes to 

(e) = -'-^J- [d^ ^ = [d^p I 

(14) 

from which one finds easily that 

(^) = T^i^gM'Y'^/XL"-'. (15) 
Collecting the results (a) to (e) we end up with (^) of the third section. 
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